function q = euler_parameter(a)
%EULER_PARAMETER Summary of this function goes here
%   Detailed explanation goes here
th = norm(a);
a = a/th;
w = cos(th/2);
xyz = sin(th/2)*a;
q = [w xyz];
end

